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ABSTRACT 



MOST time-series photometry of /j, Eri, an SB1 eclipsing binary with a rapidly- 
rotating SPB primary, is reported and analyzed. The analysis yields a number of 
sinusoidal terms, mainly due to the intrinsic variation of the primary, and the eclipse 
light-curve. New radial- velocity observations are presented and used to compute pa- 
rameters of a spectroscopic orbit. Frequency analysis of the radial- velocity residu- 
als from the spectroscopic orbital solution fails to uncover periodic variations with 
amplitudes greater than 2 kms" 1 . A Rossiter-McLaughlin anomaly is detected from 
observations covering ingress. 

^From archival photometric indices and the revised Hipparcos parallax we derive 
the primary's effective temperature, surface gravity, bolometric correction, and the 
luminosity. An analysis of a high signal-to-noise spectrogram yields the effective tem- 
perature and surface gravity in good agreement with the photometric values. From 
the same spectrogram, we determine the abundance of He, C, N, O, Ne, Mg, Al, Si, 
P, S, CI, and Fe. 

The eclipse light-curve is solved by means of EBOP. For a range of mass of 
the primary, a value of mean density, very nearly independent of assumed mass, is 
computed from the parameters of the system. Contrary to a recent report, this value is 
approximately equal to the mean density obtained from the star's effective temperature 
and luminosity. 

Despite limited frequency resolution of the MOST data, we were able to recover 
the closely-spaced SPB frequency quadruplet discovered from the ground in 2002-2004. 
The other two SPB terms seen from the ground were also recovered. Moreover, our 
analysis of the MOST data adds 15 low-amplitude SPB terms with frequencies ranging 
from 0.109 d" 1 to 2.786 d -1 . 

Key words: stars: early-type - stars: individual: ft Eridani - stars: SB1 - stars: 
eclipsing - stars: oscillations 
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1 INTRODUCTION 
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The star fj, Eri = HR 1520 = HD 30211 (B5IV, V = 4.00 
mag) w as discovered to be a single-lined spectroscopic bi- 
nary by iFrost fc Led (ll910T) at Yerkes Observa tory. An or- 
bit was derived bv lBlaauw fc van Albadal (|l963h from radial 
velocities measured on 19 spectrograms taken in November 
and December 1956 with the 82- inc h McDonald Observatory 
telescope. The orbit was revised by[mj] (|l969h who supple- 
mented the McDonald observations with the Yerkes r adial- 
velocities (published bvlFrost. Barrett fc Struvel 19*261 ) , Lick 
Observatory data of lCampbell fc Moord ( 19281 ) and Domin- 
ion As trophysical Observatory (DAO) observations of lPetrid 
(1958). Since these data span more tha n ha l f a cen tury, the 
value of the orbital period derived by iHilll (|l969T ) is given 
with six-digit precision, viz., 7.35886 d. 

The star was d iscovered to be variable in light by 
lHandler et al.l (|2004l ) while being used as a comparison star 
in the 2002-2003 multi-site photometric campaign devoted 
to v Eri. A frequency analysis of the campaign data revealed 
a dominant frequency of 0.616 d" 1 . Prewhitening with this 
frequency resulted in an amplitude spectrum with a strong 
1/f component indicating complex variation. 

Using pho t omet ric indices and the Hipparcos parallax, 
lHandler et ail (|2004T l determined the position of fi Eri in 
the HR diagram and compared it with evolutionary tracks. 
They found the star to lie close to a 6 Mq track, just be- 
fore or shortly after the end of the main-sequence phase of 
evolution. Taking into account this result, the value of the 
dominant frequency, and the fact that the u amplitude at 
this frequency is a bout a factor of two greater than the v 
and y amplitudes, lHandler et al.l (|2004r ) concluded that p 
Eri is probably an SPB star. 

[Handler et al. | (|2004 ) have also considered the possibil- 
ity that the dominant variation is caused by rotational mod- 
ulation. Assuming the dominant frequency to be equal to the 
frequency of rotation of the star and using a value of the ra- 
dius computed from the effective temperature and luminos- 
ity, they obtained ~190 kms -1 for the equatorial velocity 
of rotation, Kot- This value and the pu blished estimate s of 
Kotsini, wh ich range from 150 kms" 1 |Abt et alj|2002t ) to 
190 kms -1 H ernacca fc Perinottol 11970 )7 make rotational 
modu lation a viable hypothesis. However, lHandler et al.l 
(j2004) judged the hypothesis of rotational modulation to be 
less satisfactory than pulsations because it does not explain 
the complexity of the light variation. 

The star was again used as the comparison star for v 
Eri in the subsequent 2003-2004 multi-site photometric cam- 
paign (|jerzvkiewicz et al.ll2005l . henceforth JHS05) . The fre- 
quency of 0.616 d _1 was confirmed and five further frequen- 
cies, ranging from 0.568 to 1.206 d" 1 , were found. The mul- 
tiperiodicity, the values of the frequencies and the decrease 
of the uvy amplitudes with increasing wavelength implied 
high-radial-order g m odes, strengthening the conclusion of 
lHandler et all i|2004l ) that ^ Eri is an SPB variable. 

An attempt at determining the angular degree and 
order, I and m, of the six observed SPB term s was un- 
dertaken by iDaszvh ska-Dasz kiewicz et al.l (|2008l ). Using a 
method of mode identification devised for SPB models in 
whi ch the pulsation and rotation frequencies are compara- 
ble ( Daszvh ska-Das zkiewicz et a"l"ll2007l ; iDzie mbowsk i et al.l 
12007^ 7 these authors obtained £ = 1 and m = for the 
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dominant mode. In addition, they estimated the expected 
amplitude of this mode's radial-velocity variation to be 5.9 
kms -1 . For the remaining terms, there was some ambiguity 
in determining I and m. Nevertheless, at least four terms 
detected in the star's light variation could be attributed to 
unstable g modes of low angular degree provided that the 
star was assumed to be still in the main-sequence phase of its 
evolution. For these four terms, the predicted radial- velocity 
amplitudes ranged from 0.6 kms -1 to 13.1 kms -1 . 

In addition to confirming the SPB classification, the 
analysis of the 2003-2004 campaign data showed /i Eri to be 
an eclipsing variable. As can be seen from figs. 3 and 5 of 
JHS05, the eclipse is a probably total transit, the secondary 
is fainter than the primary by several magnitudes, and the 
system is widely detached. The eclipse ephemeris, derived 
from the combined 2002-2003 and 2003-2004 campaign data, 
can be expressed as follows: 

Min. light = HJD 2452574.04 (4) + 7.3806 (10) E. (1) 

Here, E is the number of minima after the initial epoch given 
(which is that of the middle of the first eclipse, observed in 
2002), and the numbers in parentheses are estimated stan- 
dard deviations with the leading zeroes suppressed. 

The photometric orbital period in eq uation (1) di ffers 
by more than 20 standard deviations from lHilP3 (|l969h pe- 
riod, mentioned in the first paragraph of this Introduction. 
IffiU'j (|l969h value is certainly less secure than the pho- 
tometric value because it is based on a small number of 
observations, distributed unevenly over an interval of 53.2 
ye ars. In the pe r iodo gram of these observations, computed 
by I Jerzvkiewiczl l|2009l ). the three highest peaks in the pe- 
riod range of interest, all of very nearly the same height, 
occurred at 7.35889, 7.38116 and 7 3896 2 d. The first of 
these numbers corresponds to iHuT s (1969) value, while the 
second is equal to the photometric period to within half the 
standard deviation of the latter. Since both t hese numbers 
repres ent the radial- velocity data equally well, Ijerzvkiewiczl 
(2009) concluded that the photometric period is the correct 
one. 

iBruntt fc Southworthl ()2008h observed p Eri in Febru- 
ary 2004 and August 2005 by means of the WIRE 
star tracker. After removing the SPB light- variation, they 
phased the observations wit h the orbital period of 7.381 d. 
IBruntt fc Southworthl (|2008h suggest that the eclipse depth 
is variable. 

The radial- v elocity observations of 

iBlaauw fc van Albadal dl963l ) and the photometric time- 
series data of Handler et alj l|2004h and JHS05 were used 
by Ijerzvkiewiczl (|2009l ) to revise the spectroscopic orbital 
elements, solve the eclipse light-curve, and derive the 
primary's mean density. The mean density turned out to 
be smaller than the value ob tained from the posit ion in the 
HR diagram determined by lHandler et alj (|2004f) , placing 
the star beyond the TAMS. This conflicts not only with 
the star's photometric effective temperature, but also with 
the SPB classification because pulsational instability of 
g modes disappears shortly after the minimum of T e fj in 
the m ain-sequence phase of evolution (Dzicmbowsk i et al.l 
Il993h . Improving the orbital elements of n Eri and thus 
obtaining an accurate value of the mean density of the 
primary component in order to find out how severe these 
conflicts are was the motivation for undertaking the present 



work. In addition, we expected to verify the frequencies of 
the six SPB terms found by JHS05 and discover further 
frequencies, given that the photometric amplitudes are too 
small to be detected from the ground. Finally, we planned 
to look for the SPB radial-ve l ocity variations predicted by 
Daszvhska-Daszki ewicz et alj (|2008l ) . 

In the next section, we describe the MOST times-series 
photometry of n Eri. In Section 3, we carry out a frequency 
analysis of these data and use the analysis results to re- 
move the intrinsic component of the variation, thus bringing 
out the eclipse light-curve; in this section we also improve 
the eclipse ephemeris. In Section 4, we present new radial- 
velocity observations of [s, Eri, some obtained contempora- 
neously with the MOST photometry, compute a spectro- 
scopic orbit from them, and show that in the radial-velocity 
residuals from the orbital solution there are no periodic 
terms with an amplitude exceeding 2 kms -1 . In Section 
5, we solve the eclipse light-curv e by means of the com- 
puter program EBOP (|Etze]||l98lh . Section 6 deals with the 
Rossiter-McLaughlin effect. In Section 7, we determine the 
primary component's fundamental parameters from (1) pho- 
tometric indices and the revised Hipparcos parallax and 
(2) from high-resolution spectrograms; the spectrograms are 
also used to obtain the metal abundances by means of the 
method of spectrum synthesis. In Section 8, we derive the 
primary component's mean density and discuss the star's po- 
sition in the HR diagram in relation to evolutionary tracks. 
Section 9 is devoted to assessing the SPB frequencies and 
amplitudes. The last section contains a summary of the re- 
sults and suggestions for improvements. Short-period, small- 
amplitude variations found in the MOST observations of fi 
Eri are considered in the Appendix. 



2 MOST TIME-SERIES PHOTOMETRY 

MOST (|Walker et all I2003T ) observed ^ Eri in 2009, from 
November 4 to 16. The light curve is shown in Fig. [1] The 
two wide gaps seen in the figure divide the light curve into 
three parts: (1) from the first point at HJD 2455140.1633 
to HJD 2455142.0219, (2) from HJD 2455142.5939 to 
HJD 2455146.1774, and (3) from HJD 2455146.5240 to the 
last point at HJD 2455152.1273. Part 1 throughout and part 
2 from the beginning to HJD 2455145.1394 consist of 27 
and 37 segments, respectively. The duration of most seg- 
ments is equal to about 0.045 d, but there is one shorter 
segment in part 1 (0.0265 d), and two in part 2 (0.0098 
d and 0.0356 d); the mean duration is equal to 0.0433 d. 
The mean duration of gaps between segments (not counting 
the gap between part 1 and part 2) amounts to 0.0259 d, 
and the mean interval between successive segments is equal 
to 101.4 min, the orbital period of the satellite. In part 2, 
from HJD 2455145.1554 to the end and in part 3 throughout 
there are less than one-hundred gaps, all shorter than 0.012 
d. The shortest sampling interval and the number of data- 
points per unit time differ between segments and parts. In 
all segments of part 1 and the first 28 segments of part 2, 
the shortest sampling interval is equal to 0.0004 d and there 
are about 111 data-points per hour; in the remaining nine 
segments and the non-segmented portion of part 2 and in 
part 3, the shortest sampling interval is equal to 0.0007 d 
and there are about 55 data-points per hour. 



4 M. Jerzykiewicz et al. 




1 40 1 42 1 44 146 1 48 1 50 1 52 
HJD-2455000 



Figure 1. The MOST light-curve of /i Eri. The vertical lines 
indicate estimated epochs of mid-eclipse. 
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Figure 2. The power spectrum of the out-of-eclipsc data. 

The SPB variation time-scales are much longer than 
the duration of the data segments. It is therefore surprising 
that in many segments one can see short-term variations 
with ranges reaching 5 mmag. Short-term variations can be 
also seen in the non-segmented data. Further discussion of 
these variations we defer to the Appendix. 



3 THE ECLIPSE 

3.1 Frequency analysis of the out-of-eclipse data 

In the light curve shown in Fig. [1] one can see the complex 
SPB variation and two successive eclipses. Unfortunately, in 
the first eclipse the ingress and the first half of the phase 
of totality are missing because of the first wide gap in the 
data. In addition, the eclipse is segmented (see Sect. 2), with 
the widest gap between segments (0.035 d) affecting the 
egress. We estimate the mid-eclipse epochs to be equal to 
HJD 2455142.60 ±0.03 d and HJD 2455149.98 ± 0.02 d (see 
the vertical lines in Fig.[T]). In order to obtain data suitable 
for analysing the SPB variation we modified the time series 
by omitting points falling within ±0.03 orbital phase away 
from the mid-eclipse epochs. In the process, the number of 
data points was reduced from 15687 to 14765. 

The power spectrum or periodogram of the modified 
data is shown in Fig. [2] Here and in what follows, by 'power' 
we mean 1 — a 2 (/) /a 2 , where a 2 (/) is the variance of a least- 
squares fit of a sine curve of frequency / to the data, and 
a 2 is the variance of the same data. 

As expected, the power in the periodogram in Fig. [2] is 
concentrated in the frequency range < / < 1.4 d"\ with a 
low-power alias around 14.2 d _1 , the orbital frequency of the 
satellite. This alias consists of the scaled-down pattern seen 
at low frequencies, shifted by 14.2 d" 1 , and the mirror reflec- 
tion about the / = 14.2 d _1 axis of the scaled-down pattern. 




Frequency [d ] 

Figure 3. The power spectrum of the out-of-eclipse data (top), 
after prewhitening with the frequency /i = 0.582 d _1 (middle), 
and after prewhitening with f\ and /j = 0.684 d _1 (bottom). 

The periodogram in the frequency range < / ^ 1.4 d _1 is 
displayed in the top panel of Fig. [3] The FWHM of the peri- 
odogram's highest lobe, centered at /i = 0.582 d _1 , is very 
nearly equal to 1/T = 0.083 d" 1 , where T = 12.0 d is the 
total time-span of the data. This shows that in our case the 
commonly used estimate of the frequency resolution, given 
by 1/T, should be pre ferred over the more string ent value 
of 1.5/T, advocated bv lLoumos fc Deeming! (|l978l ). 

In the middle panel of Fig. there is shown the pe- 
riodogram of the out-of-eclipse data prewhitened with the 
frequency /i. The highest lobe, centered at /2 = 0.684 d _1 , 
has already appeared as a just-resolved feature in the peri- 
odogram of the original data in the top panel of the figure. 
Finally, the bottom panel shows the periodogram obtained 
after prewhitening the data with the two frequencies, /i and 
/2. The power at low frequencies seen in this panel indicates 
a variation on a time scale longer than T. The isolated lobe 
at right is centered at 1.186 d _1 . 

Fig- H] shows the out-of-eclipse magnitudes prewhitened 
with fx, /2, and 1.186 d _1 . The long-term variation is clearly 
seen. A similar long-term variation was present in the MOST 
observations of S C eti. In the latter case it was shown by 
Ijerzvkiewic z (2007) to have been caused by inadequate com- 
pensation for background light in the MOST data reduc- 
tions, so that the maximum of the variation occurred close 
to the phase of Full Moon, and the minimum, close to the 
phase of New Moon. The present case is similar. The upper 
line in Fig. [4] represents a segment of a sine-curve having a 
period equal to a synodic month and its maximum fixed at 
the phase of Full Moon, fitted to the data together with the 
three frequencies. Also shown in the figure is a linear trend, 
fitted to the data together with the same frequencies. The 
near coincidence of the two lines proves our point. However, 
the assumption that the maximum of the long-term varia- 
tion occurs exactly at the phase of Full Moon may be unwar- 
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Table 1. Frequencies, amplitudes and phases in the MOST 
out-of-eclipse magnitudes of fi Eri. The last column contains 
the signal-to-noise ratio. 



144 146 148 

HJD-2455000 

Figure 4. The out-of-eclipse magnitudes of /i Eri prewhitened 
with /i, /2, and 1.186 d _1 . The straight line (black) is a linear 
trend and the line running above it (red) is a segment of a sine- 
curve that has a period equal to a synodic month and maximum 
at the phase of Full Moon, fitted to the data together with the 
three frequencies. (A colour version of this figure is available in 
the online journal.) 



ranted because it is not known how much of the variation 
is caused by direct illumination from the Moon (which is a 
function not only of the lunar phase, but also of the angular 
distance of the Moon from the spacecraft), and how much 
is due to increased earth-shine on moonlit nights. In addi- 
tion, assuming the variation to be strictly sinusoidal may 
be an oversimplification. Unfortunately, the total time-span 
of the data is shorter than 0.4 of a synodic month, so that 
the significance of these two factors cannot be evaluated. 
For the same reason, including the synodic frequency in the 
analysis may affect the low SPB frequencies. Therefore, we 
decided to account for the long-term variation by subtract- 
ing a linear trend. Thus, the observational equations to be 
fitted by the method of least squares to the out-of-eclipse 
magnitudes, rrn, were assumed to have the following form: 



N 

E 



A +Bti+) Aj Bia(2irfjti + $>j 



(2) 



where the second term accounts for the trend and ti = 
HJD-2455146. 

The frequencies from fz to /jv were then derived, one 
at a time, by determining the frequency of the highest peak 
in the periodogram computed after prewhitening with all 
frequencies already found and the trend. Adding the trend 
has affected the third frequency to be derived: it changed 
from 1.186 d" 1 to 1.193 d _1 . An attempt to refine the fre- 
quencies by means of the method of non-linear least squares 
was abandoned because (1) the "refined" frequencies dif- 
fered grossly from their initial values, and (2) the solu- 
tions diverged for N ^ 9. Parameters of equation (2) for 
N = 31 were computed by means of the linear least squares 
fit with frequencies read off the periodograms. The stan- 
dard deviation of the fit was equal to 0.64 mmag, the first 
term, Aq = 4297.843 ± 0.008 mmag, and the second term, 
B = 0.6253 ± 0.0026 mmagd -1 . The frequencies, fj, the am- 
plitudes, Aj, and the phases, are listed in Table 1 in or- 
der of decreasing amplitude. The signal-to-noise ratios, S/N, 
given in the last column, were computed using Aj as S, and 
as TV, the mean amplitude over a 1 d -1 interval around fj in 
the amplitude spectrum of the data prewhitened with all 31 
frequencies. The decrease of S/N with Aj is not monotonic 
because TV is a function of fj. The amplitudes, and there- 
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fore the order of terms in Table 1, depend on the number of 
terms taken into account. This is why the third frequency 
derived from the periodograms is fi in the table. The fre- 
quencies which differ by less than 1/T from a frequency of 
larger amplitude are indicated with an asterisk. We included 
these frequencies because omitting them would defy our pur- 
pose of removing as much out-of-eclipse variation as possible 
from the light-curve. By the same token, we included terms 
with S/N < 4.0 and those with fj > 15 d^ 1 although the 
former are in s ignific ant according to the popular criterion of 
iBreger et alj (| 19931 ) while the latter are probably artefacts 
because they have frequencies that are either close to the 
whole multiples of the orbital period of the satellite or dif- 
fer from the whole multiples by ~1 d _ . As shown in the 
Appendix, these frequencies are related to the short-term 
variations mentioned in Sect. 2. 



3.2 The eclipse light-curve and the corrected 
ephemeris 

The residuals, O — C, where O is the observed magnitude 
and C is the magnitude computed from a least-squares fit 
of equation (2) with TV = 31 to the out-of-eclipse data, are 
plotted in Fig. [5] From the residuals, the mid-eclipse epoch 
of the second eclipse could be derived much more accurately 



6 M. Jerzykiewicz et al. 




1 40 1 42 1 44 146 1 48 1 50 1 52 
HJD-2455000 



Figure 5. The residuals, O — C, from a least-squares fit of equa- 
tion (2) with N = 31 to the out-of-eclipse data. 

than was possible with the original data (see the first para- 
graph of Sect. 3.1). The result is HJD 2455149.973 ± 0.002. 
Using this number and equation (1) we arrive at the follow- 
ing corrected ephemeris: 

Min. light = HJD 2455149.973 (2) + 7.38090 (12) E, (3) 

where the notation is the same as before, except that the 
minima are now counted from the new, more accurate initial 
epoch. Note that the new value of the orbital period differs 
from the old value by about 0.3 of the standard deviation of 
the latter. 



4 THE SPECTROSCOPIC ORBIT 
4.1 The radial- velocity data 

On 24 nights between 2 October 2009 and 8 February 2010, 
HL, MHa, and MHr obtained 105 spectrograms of \i Eri 
with the coude echelle spectrograph at the 2-m Alfred Jen- 
sch telescope of the Thiiringer Landessternwarte Tautenburg 
(TLS). The spectral resolving power was 63,000. The spec- 
trograms were reduced using standard MIDAS packages. Re- 
ductions included bias and stray light subtraction, filtering 
of cosmic rays, flat fielding, optimum extraction of echelle 
orders, wavelength calibration using a ThAr lamp, correc- 
tion for instrumental shifts using a large number of O2 tel- 
luric lines, normalization to local continuum, and weighted 
merging of orders. 

On eight nights, 4 to 14 November 2009, JM-Z obtained 
43 spectrograms of /i Eri with the 91-cm telescope and the 
fiber-fed echelle spectrograph FRESCO of the M.G. Fra- 
castoro station of the Osservatorio Astrofisico di Catania 
(OAC). The spectral resolving power was 21,000. The spec- 
trograms were reduced with 

Finally, on eight nights between 29 September 2009 and 
22 January 2010, WD obtained 64 spectrograms of the star 
with the Poznari Spectroscopic Telescopqj (PST). The re- 
ductions were carried out by MF and KK by means of the 
Spectroscopy Reduction Packag^f). 

Radial velocities were determined from the spectro- 
grams by EN by means of cross-correlation with a template 

1 IRAF is distributed by the National Optical Astronomy Obser- 
vatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc. 

2 www.astro.amu.edu.pl/PST/ 

3 http://www. astro. amu.edu. pl/~chrisk/MOJA/pmwiki.php?n= 
SRRSRP 



Table 2. Nightly mean radial- velocities of fi Eri. 



(HJD) 


(RV> 
[kms -1 ] 


(S/N) 


N 


Obs 


94^1 04 fi4Q1 


28 35 


87.7 


3 


PST 


24551 06 591 3 


—8.02 


191.5 


2 


TLS 


94^1 07 fi1 Q4 


13.82 


179.4 


■ 


TLS 


94^51 1 3 fi1 SO 


— 7.24 


210.0 


I 


TLS 




—4.29 


200.5 


2 


TLS 


94^51 ^(S 691 


—4.76 


210.0 


1 


TLS 




36.22 


91.7 


4 


OAC 


2455141 47fi1 


30.84 


110.5 


4 


PST 


94SS141 Sfifil 


30 12 


226.0 


2 


TLS 


2455142 5027 


13.82 


180.0 


1 


TLS 


94SS1 43 S04& 


— 11.36 


150.0 


1 


TLS 


24 c i c S1 5660 


— 10.49 


94.8 


5 


OAC 


94SS1 44 4SS r i 


4.17 


110.0 


I 


OAC 


24551 4fi 51 52 


32.59 


87.6 


2 


OAC 


94SS1 47 SS7& 


36.65 


216.0 


2 


TLS 


94^5147 57fi7 


37.65 


96.4 


7 


OAC 


94SS1 48 S4&4 


35 66 


92.5 


7 


OAC 


2455149 5758 


21.62 


88.3 


g 


OAC 


24551 50 5478 


—5.55 


100.7 


9 


OAC 


24551 =i1 4388 


—2.40 


75.3 


4 


PST 


04^ "SI 4Q78 


37.30 


72.6 


21 


PST 


24551 55 5034 


37.10 


215.1 


2 


TLS 


24551 56 51 90 


30.64 


67.3 


g 


PST 


24551 57 41 54 


9.53 


69.0 


3 


PST 


24551 58 5251 


—9.05 


240.0 


1 


TLS 


2455161 4549 


37.23 


87.7 


14 


PST 


2455162.4466 


32.40 


180.0 


1 


TLS 


2455163.4780 


32.42 


170.0 


1 


TLS 


2455168.4333 


31.12 


190.0 


1 


TLS 


2455170.4825 


37.45 


255.1 


2 


TLS 


2455175.3956 


29.30 


230.0 


1 


TLS 


2455192.4014 


36.54 


215.1 


4 


TLS 


2455194.3866 


2.94 


230.0 


1 


TLS 


2455199.3759 


37.07 


179.8 


15 


TLS 


2455201.4046 


21.68 


225.1 


4 


TLS 


2455202.3392 


-12.10 


191.4 


34 


TLS 


2455219.3368 


25.07 


76.0 


7 


PST 


2455228.3739 


39.52 


260.0 


1 


TLS 


2455231.2637 


3.63 


270.0 


1 


TLS 


2455236.2367 


34.70 


200.0 


1 


TLS 



spectrum in the logarithm of wavelength. A synthetic spec- 
trum, computed assuming the photometric T e s and logg 
(see Sect. 7.1), a microturbulence velocity of 4.0 kms -1 , 
Kotsini = 136 kms -1 , solar helium abundance, and the 
EN metal abundances from Table 6, served as the template. 
The wavelength intervals used in the cross-correlations in- 
cluded He I lines and lines of the elements listed in Table 6, 
had high signal-to-noise ratio (S/N), and were not affected 
by telluric lines. The S/N were determined with the IRAF 
task splot (function "m") from the continuum in the 5100 
A to 5800 A range which contains only very faint lines. The 
S/N (per pixel) of the spectrograms ranged from 70 to 270 
for TLS, from 60 to 110 for OAC, and from 40 to 120 for 
PST; the overall mean S/N values were 198 ± 3, 93 ± 2, and 
78 ±2 for TLS, OAC, and PST, respectively. 
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Table 3. Orbital elements of fi Eri computed from the radial- 
velocities of Table 2, with the systematic corrections added to 
the OAC and PST data (see text and Fig.0. 



P olh (assumed) 7.38090 d 

T HJD 2455143.254 ± 0.067 d 

u> 160.5° ±4.5° 

e 0.344 ±0.021 

7 20.58 ±0.43 kms" 1 

Ki 24.24 ±0.53 km s" 1 

aisini 3.320 ± 0.077 R Q 

f(M) 0.00902 ± 0.00063 M 
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Figure 6. The TLS, OAC, and PST radial-velocities from Ta- 
ble 2, with the systematic corrections in the cases OAC and PST 
as indicated, plotted as a function of orbital phase. Also shown is 
the radial- velocity curve (line), computed from the spectroscopic 
elements of Table 3. 

4.2 The spectroscopic elements 

The number of spectrograms taken on a night ranged from 
one to 34. Using the individual radial-velocities of the pre- 
ceding section for deriving a spectroscopic orbit would give 
excessive weight to the few nights with a large number of 
observations. To avoid this, we took weighted means of the 
individual velocities on nights with two or more observations 
using the S /N of the individual spectrograms as weights; the 
epochs of observations and S/N were likewise averaged. In 
Table 2 are given the single (N = 1) or the nightly weighted 
mean (N ^ 2) epochs of observations, (HJD) , the radial ve- 
locities, (RV), and the signal-to- noise ratios, (S/N). In the 
case of JD 2455201, the means were computed from the four 
points unaffected by the Rossiter- McLaughlin effect (see be- 
low, Sect. 6). Using the (HJD), (RV), and (S/N) from Ta- 
ble 2, we computed a spectroscopi c orbit by means o f the 
non-linear least squares method of Schlesi nger (1908|). We 
applied weights proportional to (S/N). The orbital period 
was assumed to be equal to the photometric period given 
in equation (3). An examination of the residuals from the 
solution revealed small systematic shifts between the TLS, 
OAC, and PST radial- velocities. Corrections to the OAC 
and PST radial- velocities which minimized the standard de- 
viation of the solution turned out to be —0.09 ± 0.73 km s -1 
and -0.15 ±0.71 kms -1 for OAC and PST, respectively. 
The elements of the orbit, computed with these corrections 
taken into account, are listed in Table 3. Figure [6] shows 
the TLS radial-velocities and the corrected OAC and PST 
radial- velocities, plotted as a function of phase of the orbital 
period. The solid line is the radial-velocity curve computed 
from the elements listed in Table 3. 



4.3 Looking for intrinsic radial-velocity variations 

In an attempt to detect intrinsic radial-velocity variations, 
we examined the residuals from the orbital velocity curve 
of the 40 nightly mean radial- velocities used to derive it. 
A periodogram of these data, computed with weights pro- 
portional to (S/N) (see Table 2), showed the highest peak 
at / = 0.812 d - . Interestingly, this frequency is close 
to the JHS05 frequency (/a in th eir notation) for which 
Daszvhsk a-Daszkiewicz et al.l l|2008l ) predict the largest ve- 
locity amplitude (see the Introduction). However, the pre- 
dicted amplitude of fi, was equal to 13.1 kms -1 , while the 
amplitude of / is equal to only 1.9 ±0.4 kms -1 . In addi- 
tion, a periodogram of the residuals from the orbital ve- 
locity curve of the individual radial-velocity measurements 
(without those affected by the Rossiter-McLaughlin effect) 
showed the highest peak with the same amplitude but at a 
frequency of 0.626 d _1 , unrelated to /. We conclude that 
both these frequencies, / and 0.626 d _1 , are spurious, and 
that in the radial-velocity residuals there are no periodic 
terms with an amplitude exceeding 2 kms -1 . 

Unlike the PST and TLS observations, which fall mainly 
outside the interval covered by the MOST observations 
of n Eri, all 43 OAC radial-velocities were obtained con- 
temporaneously with the MOST photometry. Taking ad- 
vantage of this, we have looked for a correlation between 
the OAC radial-velocity residuals and the trend-corrected 
MOST magnitudes. No correlation was found. Thus, the 
largest deviations from the orbital velocity curve, amounting 
to about ±4.7 kms -1 (see Fig. [6}, are probably instrumen- 
tal. 



5 THE ECLIPSE LIGHT-CURVE SOLUTION 

For computing the p aramete rs of [i Eri from the eclipse 
light-curve we used lEtzel'd i| 19811) computer program 
EBOP, which is based on the Nelson-Dav is-Etzel model 
|Nelson fc Davislll972l ; |popper fc Etzellll98 ll ). The program 
is well suited for dealing with detached systems such as 
the present one. The spectroscopic parameters uj and e, 
needed to run the program, were taken from Table 3. The 
primary's limb-darkening coefficient we com puted by con- 
volvi ng the MOST KG filter transmission dWalker et al.l 
2003) with the monochromatic limb-darkening coefficients 
for T efI = 15 670 K and log 3 = 3.5 (see Section 7.1) 
and [M/H] = (Section 7.2). The monochromatic limb- 
darkening coefficients were interpolated in the tables of Wal- 
ter Van Hammfl For the primary, the result was 0.308. 
For the secondary, assuming the effective temperature to 
scale as the square root of the radius and taking the ratio 
of the radii to be equal to 0.135 (see Table 4), we obtained 
Tcft = 5 758 K. Using this value, we found the central sur- 
face brightness of the secondary in units of that of the pri- 
mary (which EBOP uses as a fundamental parameter) to 
be equal to 0.018. Assuming further the solar value of 4.4 
for logg and solar metallicity, and using the same tables 
of monochromatic limb-darkening coefficients as above, we 
computed the secondary's limb-darkening coefficient to be 

4 http: //www2. fiu .edu/~vanhamme/limdark. htm, see also 
IVan Hammel jl993h 
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Figure 7. Normal points computed in adjacent intervals of 0.0005 
orbital phase using the O — C residuals shown in Fig. [5] The solid 
line is the synthetic light-curve, computed with the parameters 
(4) of Table 4. 



Table 4. The parameters of fj, Eri derived by means of EBOP. Data 
(1) ... (4) are explained in the text. SE is the EBOB standard error 
of one observation in mmags. 



Data 


SE 


i 


k 


r P 


(1) 


0.71 


82.06±0.07 


0.13507±0.00011 


0.2094±0.0005 


(2) 


0.48 


82.11±0.12 


0.13524±0.00021 


0.2092±0.0008 


(3) 


0.64 


82.42±0.07 


0.13502±0.00010 


0.2072±0.0004 


(4) 


0.41 


82.44±0.11 


0.13539±0.00018 


0.2071±0.0007 



equal to 0.628. The remaining parameters we set to their 
default EBOP values. The unknowns to be solved for were 
the orbital inclination, i, the ratio of the radii, k = R s /R p , 
and the radius of the primary, r p = R p /a, in units of the 
semi-major axis of the relative orbit, a. Initial values of the 
unknowns were derived by trial and error. Tests showed that 
the results were insensitive to the initial values. 

The results depend slightly on the light-curve data used. 
In Table 4 are shown the parameters obtained using (1) the 
O — C residuals from a least-squares fit of equation (2) with 
N = 20 to the out-of-eclipse data, i.e., without terms fj > 
15 d _1 , (2) normal points computed from the O—C residuals 
in adjacent intervals of 0.0005 orbital phase; (3) and (4): 
the same as (1) and (2), respectively, but with N = 31. A 
synthetic light-curve computed with the parameters derived 
from data (4) (see Table 4) is compared with the data in 
Fig. [7] for orbital phases from —0.1 to 0.1, the synthetic 
light-curve is compared with the data in Fig. [8] 

As can be seen from Fig. the fit of the computed 
light-curve (line) to the observations (points) is not perfect 
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Figure 8. The same as Fig. [7] but for a limited interval of orbital 
phase. 
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Figure 9. Upper panel: individual radial-velocity measurements 
on JD 2455201 (open circles) plotted as a function of orbital 
phase. The short vertical lines indicate the epochs of first and 
second contact. Also shown are the radial-velocity curve, com- 
puted from the spectroscopic elements of Table 3 (solid line), 
and the radial-velocity curve plus the anomaly due to the 
Rossiter-McLaughlin effect, computed as explained in the text for 
Vrot sini =115 kms" 1 (dashed line). Lower panel: radial- velocity 
residuals (open circles) and the Rossiter-McLaughlin anomaly 
for Fj-otsini = 115 kms" 1 , shifted upwards by 4.67 kms -1 
(dashed line), bracketed by those for V ro t sini = 70 kms -1 and 
Vrot sini = 160 kms -1 (the lower and upper dotted line, respec- 
tively) . 



at the bottom of the eclipse: before mid-eclipse the points 
fall below the lines, i.e., the system is fainter than computed, 
while after mid-eclipse the reverse is true. This looks like a 
signature of the rotational Dopp ler beaming, i.e., photomet- 
ric R ossiter-McLaughlin effect (jGroolJ 120121: IShporer et alJ 
|2012| ). However, from equation (3) of lShporer et alJ ( 20121 ) 
it follows that in the present case the amplitude of the effect 
would be an order of magnitude smaller than the observed 
deviation. 

The depth of the secondary eclipse predicted by the 
EBOP solution is 0.3 mmag, much smaller than the noise 
seen in Fig. [7] In the infrared, the depth of the secondary 
eclipse will be greater because the contribution of the pri- 
mary to the total light of the system would decrease with 
increasing wavelength. The predicted MOST magnitude dif- 
ference between components is 8.8 mag. Assuming that this 
number represents the magnitude difference in the V band 
and using the standard V — J and V — K colour indices, one 
gets 2 mmag and 3 mmag as the depths of the secondary 
eclipse in J and K, respectively. Thus, detecting the sec- 
ondary eclipse of fj, Eri would be a task for a future infrared 
satellite. 



6 THE ROSSITER-MCLAUGHLIN EFFECT 

On JD 2455201, there were 25 measurements, the first one at 
HJD 2455201.3958, and the last at HJD 2455201.5034. Ac- 
cording to equation (3), the corresponding orbital phases 
(reckoned from the epoch of mid-eclipse) are equal to —0.033 
and —0.018. Thus, the measurements start shortly before 
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the ingress begins, and last to just after it ends. This is il- 
lustrated in the upper panel of Fig. [9] where the epochs of 
first and second contact, predicted by the eclipse solution 
(4) of Table 4, are indicated by short vertical lines. Also 
shown in the figure is a fragment of the orbital velocity curve, 
computed from the spectroscopic elements of Table 3 (solid 
line in the upper panel), and the deviations from this curve 
due to the Rossiter-McLaughlin effect (dashed line), com- 
puted assuming (1) the direction of rotation concordant with 
the direction of orbital motion, and (2) Kot sin i of the pri- 
mary component equal to 115 kms -1 , a value that gave the 
best least-squares fit to the data. The lower panel of Fig. [9] 
shows the fit to the residuals, bracketed by the Rossiter- 
McLaughlin deviations for Kot sini equal to 70 kms -1 and 
160 kms -1 (the lower and upper dotted line, respectively). 
In computing the deviations we used formulae developed by 
IPetrid (|l938l ). These formulae assume the axis of rotation 
of the star to be perpendicular to the orbital plane, neglect 
limb darkening, and are valid for circular orbits. The last 
assumption was moderated to first ord er in th e eccen tricity 
by means of the prescription given by iRusselll l|l912h . Tak- 
ing into account the effect of limb darkening w ould reduce 
the computed deviations (see lOhta et al.|[2005l ). so that the 
best fitting Kot sin i would be greater by a factor of ~1.5. 
We conclude that the increase of the radial- velocity residuals 
from about 4.5 kms -1 to about 6 kms -1 , seen in the lower 
panel of Fig. is consistent with Kot sin i of about 170 ± 45 
kms -1 , but we postpone a more detailed discussion of the 
Rossiter-McLaughlin effect in /i Eri until more accurate data 
become available. Note, however, that although the value of 
Kot sin i is not well constrained at this stage, the direction 
of rotation is determined beyond doubt. 

On JD 2455201, the mean residual from the orbital 
solution amounted to 5.49±0.15 kms -1 . Subtracting the 
Rossiter-McLaughlin deviations (dashed line in Fig. [9]) re- 
duces this value to 4.67 ± 0.10 kms -1 . The residual is prob- 
ably instrumental (see Sect. 4.3). 



7 STELLAR PARAMETERS 

7.1 From photometric indices and the revised 
Hipparcos parallax 

The colour excess of fi Eri is small but not negli gible. Using 
the st ar's B — V and U — B colour indices from iMermilliodl 
|l99ll ) and lJohnson'sl l|l966f ) standard two-colour relation for 
main- sequence stars, we get E(B — V) = 0.018 ± 0.006 mag. 
From the Stromgren uvby indices (|Hauck fc Mermilliodl 
1 19981 ) and ICrawford'sl l| 19781 ) (b - y) - c relation, we 
have E(b - y) = 0.013 ±0.004 mag or E(B - V) = 
0.018 ± 0.005 mag. The der eddened colours and the /3 index 
l|Hauck fc Merm illiod 1998) then yield the effective temper- 
ature, T e fr, the bolometric correction, BC, and the logarith- 
mic surface gravity, \ogg. The results are listed in Table 5, 
in columns three, four, and five, respec tively. The BC val- 
ues obtained using the calibrations of iBalonal (1 19941 ) and 
iFlowerl (| 19961 ) were adjusted slightly (by 0.03 mag and 0.01 
mag, respectively) in order to make them con sistent with 
BCf^ = —0.07 mag, the zero point ad opted by lCode et al.l 
l| 19761 ). As can be seen from Table 5, iFlower'sl <|l996D T eS 
and BC are much lower than the other ones. The for- 



Table 5. Fundamental parameters of fi Eri obtained from the photo- 
metric indices listed in column one, using the calibration referenced 
in column two. 



Data 


Calibration 


T e ff 


BC 


log 9 


(B - V) 


Code et al. (1976) 


15 700 


-1.47 




(B - V) 


Flower (1996) 


14 610 


-1.17 




CO 


Davis k. Shobbrook f 1977) 


15 670 


-1.46 




CO, 


UVBYBETA^ 


15740 




3.44 


co, |8 


Balona (1994) 


15 580 


-1.28 


3.52 



a A FORTRAN progra m based on the grid published by 
iMoon fc Dworetskvl l ll985h . Written in 1985 by T.T. Moon of the Uni- 
versity London and modified in 1992 and 1997 by R. Napiwotzki of 
Universitaet Kiel. 



mer values were obtained by interpolation from t his au- 
thor's table 3. Using the polynomial fits provided bv lFlowerl 
l|l996h and corrected by iTorres! i|2010T ) results in T eS and 
BC virtually identi cal with those given in Table 5. Omit- 
ting IFlower'sl (Il996l ) values, we get 15 670 K and -1.40 mag 
for the mean T c h and BC, resp ectively. Our value o f T e fr 
is identical to those derived by lHandler et al"] (|2004l ) and 
iDaszvriska-Daszkiewicz et al. ll200St) fro m Geneva photome- 
trv. However. while lHandler et aiT (2004) estimated the stan- 
dard deviation of their value to be equal to 100 K, we believe 
that the standard deviation of T G fi cannot be smaller than 
the unc ertainty of temperature calibrations, i.e., about 3 
percent (|Napiwotzki et al.lll993l;|Jerzvkiewicdll994l) or 470 K 
for the T c ft in question. IPaszvhska-Daszkiewic Terall B 
give a somewhat more optimistic value of 365 K. 

Taking into account the uncertainty of calibrations, we 
estimate the standard deviation of the bolometric correction 
to be equal to 0.08 mag. Finally, the mean value of logg is 
equal to 3.5 dex (see Table 5), with an estimated standard 

deviation of 0.15 dex. 

;From the revised Hip parcos paral l ax |van Leeuwenl 

I2007I ). the V magnitude from IMermilliodl i|l99ll ). E(B - V) 



from the first paragraph of this section, and assuming Ry = 
3.1, we get M v = -2.06 ±0.07 mag, M bo i = -3.46 ±0.10 
mag, and logL/LQ = 3.280 ±0.040. The standard devia- 
tions of M b oi and log L/Lq are determined by the standard 
deviation of BC and the standard deviation of My (in that 
order). In computing logL/LQ, we assumed M bol Q = 4.74 
mag, a value consistent with BCf^ = —0.07 mag we adopted 
and Vq = -26.76 mag (lTorresll201ol ). 



7.2 From TLS spectra 

A Fo urier-transform based program KOREL (|Hadraval 
20063, b) was used to decompose the stellar spectrum and 
the telluric contributions with the orbital elements obtained 
in Section 4.2 adopted as starting values and the orbital 
period fixed at the photometric value. The result was a 
high-S/N spectrum of fj, Eri, free of telluric contributions, 
representing a mean spectrum, averaged over all RV-shifted 
spectra of the TLS time-series. The analysis of this mean 
spectrum was carried out in dependently by HL and EN. HL 
used the program SynthV l|Tsvmball (l996l ) to compute a 
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grid of synthetic spectra from LTE m odel atmospheres ca l- 
culated with the program LLmodels JShulyak et al.ll2004) 
For a detailed description of the method see lLehmann et all 
(2011). For the analysis, the full wavelength range 4720 A to 
7300 A was used. The atomic line parameters were adopted 
from the Vienna Atomic Line Database (VALD)0- HL has 
derived the following parameters: T e g = 15 590 ± 120 
K, log g = 3.55 ±0.04, the microturbulence velocity, £ = 
0.0 ±0.4 kms -1 , and Kotsini = 130 ±3 kms -1 . For the 
He abundance, HL has obtained 1.3 ±0.1 times the solar 
va lue. The errors a r e la e rrors from \ 2 statistics as described 
in lLehmann et all (|201lh . They represent a measure of the 
quality of the fit but do not account for imperfections of 
the model. Small corrections to the local continua of the ob- 
served spectra have been applied by comparing them to the 
synthetic spectra in the process of parameter determination. 
Possible correlations between the parameters are included by 
deriving the optimum values and their errors from the mul- 
tidimensional surface in x 2 of the correlated parameters. A 
careful check for ambiguities between the parameters was 
performed. Besides the well-known correlation between T e s 
and [M/H], we observe in the given temperature range a 
slight degeneracy between T e s and log g and a strong corre- 
lation between logy and the He content. For the latter, we 
derive dlogg/dHe = —0.15 from our LLmodel atmospheres, 
where He is the He abundance in units of the solar value. 
An enhanced He abundance generat es additional pressure 
jAuer et al.lll966l : lFossati et al.ll2009l ). changing in our case 
log g from 3.36 for solar He abundance to 3.55 for 1.3 times 
sol ar abundance . The e ffect on T e ff is marginal, as also found 
by iFossati et al.l (|2009l ). We determine dT eff /dHe = 22 K. 
No degeneracy between any of the other parameters was 
found. A test of the influence of the abundances of different 
chemical elements showed that only a modification of the 
He, Ne, Si, S, and Fe abundances had a significant effect on 
the quality of the fit. Except for Ne, the metal abundances 
turned out to be very nearly solar. The numbers are given in 
the second column of Table 6; the standard deviations were 
estimated from \ 2 statistics as described before. 

EN used the same mean spectrum as HL did, a grid of 
line-blanketed LTE model atm ospheres computed with the 
ATLAS9 code l|Kurucdll993al l and synthetic spectra com- 
puted with the SYNTHE code (|Kurucdll99"3bh . The stellar 
line identification in the e ntire observed spectral ra nge was 
done using the line list of ICastelli fe Hubrid (J2004J). To be- 
gin with, the photometric T e g and log g (see Sect. 7.1) were 
verified by comparing the observed profiles of the H,9 and 
Ha lines with the synthetic profiles, computed for T e fr from 
15 000 K to 16 000 K with a step of 50 K, and log g from 3.0 
dex to 4.0 dex with a step of 0.1 dex. The best match was ob- 
tained for 15500 K < T cft < 15700 K and 3.5 dex < logg < 
3.6 dex, in excellent agreement with the photometric val- 
ues and with HL. Then, parts of the spectrum suitable for 
the analysis were selected. Next, the elements which influ- 
ence every chosen part were identified. Finally, assuming 
the photometric T e g and logg, £ = 4.0 kms -1 , and solar he- 
lium abundance, the values of V TO t sin i and the abundances 
were determined by means of the spectrum synthesis method 



5 http : / /ams. astro. univie.ac. at / ^vald, see also iKupka et al] 




4840 4860 4880 4900 

Wavelength [A] 




6480 6523 6565 6608 6650 

Wavelength [A] 

Figure 10. Upper panel: the observed mean profile of the H/3 
line (black, solid) and the synthetic profile, computed for T e ff = 
15 670 K, logg = 3.5 dex, £ = 4.0 kms" 1 , and Viotsini = 
136 kms -1 (red, dashed). Lower panel: the same for the Ho line. 
(A colour version of this figure is available in the online journal.) 
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described in iNiemczura et al.l (|2009l ) . The rotational veloc- 
ity was derived from every part of the spectrum, whereas 
the abundance of a given element was determined from all 
parts in which this element was identified. The final values 
of V^-otsini, and of the abundances are average values ob- 
tained from all parts in which they were determined. The 
final value of Kotsini is equal to 136 ±6 kms -1 . A com- 
parison of the TLS mean profiles of the H/3 and Ha lines 
with the synthetic profiles, computed with the photometric 
T e ff and logp, £ = 4.0 kms -1 , and Kotsini = 136 kms -1 , 
is presented in Fig. [10] 

The abundances with their standard deviations are 
given in the third column of Table 6, and the number of 
parts of the spectrum from which they were derived, n, in 
column four. T he solar abundances li sted in the last column 
are quoted from lAsplund et all (|2006l l. Both, HL and EN de- 
termined the abundances of Ne, Si, S and Fe. Except for Si, 
the agreement between the HL and EN abundances is within 
la. In the case of Si, HL's logiV is virtually identical with 
log Nq. As can be seen from Table 6, the C, N and O abun- 
dances determined by EN are slightly smaller then solar, 
whereas Al (only one part of the spectrum used) and P are 
enhanced. The abundances of Mg, S, CI, Fe and Ni obtained 
by EN are close to the solar values. On the other hand, 
Ne is overabundant by 0.71 ± 0.11 dex according to EN, and 
by 0.56 ±0.15 dex according to H L. The Ne ab u ndanc es of 
B-type stars were inv estigated by ICunha et all l|2006h and 
iMorel fe Butlerl (|2008l l. In these papers a non-LTE approach 
was adopted; the Ne abundances turned out to be equal 
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Table 6. The elemental abundances of /lEri derived by HL 
and EN; n is explained in the text. 



Element 


log 


\ogN 


n 


log7V 




F 1 1 






He 


11.04±0.034 






10.93 


C 




8.26±0.27 


6 


8.39 


N 




7.41±0.29 


3 


7.78 







8.49±0.31 


3 


8.66 


Ne 


8.40±0.26 


8.55±0.11 


13 


7.84 


Mg 




7.51 


2 


7.53 


Al 




6.77 


1 


6.37 


Si 


7.49±0.21 


7.11±0.12 


7 


7.51 


P 




5.82±0.21 


5 


5.36 


S 


6.94±0.20 


7.06±0.18 


18 


7.14 


CI 




5.51 


1 


5.50 


Fe 


7.45±0.30 


7.34±0.23 


20 


7.45 


Ni 




6.21±0.15 


3 


6.23 



to 8.11 ±0.04 and 7 .97 ± 0.07 dex, respectively. In addition, 
ICunha et al] (2006) found that the difference between the 
LTE and non-LTE abundance derived from the Ne I 6402 A 
line amounted to 0.56 dex. This result is close to the differ- 
ence between the HL and EN LTE-de termi nations and the 
non-L TE values of ICunha etafl (|2006h and lMorel fe Butlerl 
{200$) . 

In addition to differences in methodology and software, 
an important difference between the HL and EN analyses 
is the value of the microturbulence velocity. While HL de- 
rived £ = O.Okms -1 , EN assumed 4.0 kms -1 following 
iLefever et all fcOlOft . The main consequence of this differ- 
ence is the following: For a strong line, greater £ results 
in a greater equivalent width of the line in the synthetic 
spectrum, so that matching the synthetic and observed line- 
profiles requires a smaller abundance than would be the case 
for a smaller £. The determination of the Si abundance is 
dominated by the strong Sill doublet at 6347/6371 A. This 
explains why EN obtained appreciably smaller abundance of 
Si than did HL (see Table 6). In the case of Fe, there was one 
strong line, viz., Fel at 5057 A, and dozens of fainter lines; 
consequently, the effect is smaller than for Si. The effect is 
absent in the case of Ne and S because no strong lines were 
used in the analysis. 



8 MEAN DENSITY AND THE HR DIAGRAM 

Using the effective temperature and luminosity obtained in 
Sect. 7.1 and assuming a mass of 6.2 ±0.2 Mq (see below) 
we find the mean density of the primary component, (p) p h = 
0.0418 ± 0.0057 gem -3 , where the suffix "ph" is a reminder 
of the photometric provenance of this value. Surprisingly, the 
largest contribution to the large standard deviation does not 
come from the standard deviation of T e g, but from that of 
L: if the latter were omitted, the standard deviation of (p) p h 
would be reduced to 0.0033 gem -3 ; omitting the former 
yields a standard deviation of 0.0048 gem -3 . 

A value of the mean density independent of photometric 
calibrations can be derived in the present case because for 
a binary with the primary's mass much larger than that of 



the secondary, the spectroscopic and photometric elements 
can be combined to derive a mass-radius relation 

M p ~ R 3 p +X , (4) 

where R p is the radius of the primary and i is a small 
positive number, slowl y varying with M p , the assu med 
mass of the primary |Pigulski fe Jerzvkiewiczl Il98ll ). A 
consequence of this relation is that the mean density of 
the primary computed from the orbital solution, (p) or b, 
is weakly dependent on M p . This fact has been used by 
iDziembowski fc Jerzvkiewiczl (1 19961 ) to derive a mean den- 
sity of the p Cephei primary of 16 Lac (EN Lac), an SB1 
and EA system rather similar to p Eri. In the present case, 
the spectroscopic parameters K\ and e from Table 3 and 
the photometric parameters r p and i from line (3) or (4) 
of Table 4 give (p) orb = 0.0345 gem -3 , 0.0348 gem -3 , and 
0.0350 gem -3 for M p = 5 Mq, 6 Mq, and 7 Mq, respec- 
tively, with the standard deviations equal to 0.0032 gem -3 . 
The last value is determined mainly by the standard devi- 
ation of K\ : if the standard deviation of K\ were omitted, 
the standard deviation of (p) or b would be equal to 0.0012 
gem -3 ; the contributions of the standard deviations of the 
photometric parameters are negligible. For M p = 6.2 Mq, 
the difference between the two values of mean density, 
A (p) = (P)ph - (p)orb, is equal to 0.0070 ± 0.0065 gem -3 . 
If the photometric parameters from line (1) or (2) were 
used, (p)orb would be 0.0010 gem -3 smaller than those given 
above and A(p> would amount to 0.0080 ± 0.0065 gem -3 . 
In what follows we use the larger values of (p)orb, i-e., those 
obtained with the photometric parameters from line (3) or 
(4) of Table 4. 

Figure [TJJ shows an HR diagram in which p Eri is plot- 
ted using the effective temperature and luminosity obtained 
in Section 7.1. Also shown are X = 0.7, Z — 0.02 evolution- 
ary tracks for 5 Mq, 6 Mq, and 7 Mq. The tracks were 
kindly computed for us by Dr. J. Daszyriska-Daszkiewicz 
with the OPAL opacities, A04 mixture, no convective-core 
overshooting, and V TO t = 160 k ms -1 on the ZAMS. The code 
was th e same as that used bv lDaszvriska-Daszkiewicz et al.l 
(2008). The position of the star in relation to the evolu- 
tio nary tracks is very n early the same as that in figure 5 
of lHandler et al l (|2004l ). Thus, these authors' conclusion 
that p Eri is just before or shortly after the end of the 
main-sequence phase of evolution (see the Introduction) re- 
mains unaltered. As can be seen from Fig. 1111 the assump- 
tion that the star is still on the main sequence leads to 
M p = 6.2 ±0.2 M . 

The thick solid line in Fig. [TT] connects points for which 
(p)orb is equal to the evolutionary tracks' mean density. The 
dashed lines run at (p) or b±2a. For M p > 5 Mq, the line 
crosses the tracks at an advanced stage of hydrogen burning 
in the core, a phase of secondary contraction, or a phase 
of shel l hydrogen-burning. This is the same conclusion as 
that of lHandler et al.l (|2004l ) but independent of photometric 
calibrations and the Hipparcos parallax. 

9 THE SPB TERMS 

9.1 A comparison with the ground-based results 

The six SPB terms, derived from the 2002-2003 and 2003- 
2004 multi-site photometric campaigns (MSC, see the In- 
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Figure 11. fi Eri in the HR diagram. The point with error bars 
was plotted using the photometric effective temperature and the 
luminosity computed from the revised Hipparcos parallax and the 
photometric BC (see Section 7.1). The thin solid lines are X = 
0.7, Z = 0.02 evolutionary tracks for 5 Mq, 6 Mq, and 7 Mq. The 
thick solid line connects points for which the primary component's 
mean density derived from the orbital solution, (p) rb> is equal 
to the evolutionary tracks' mean density. The dashed lines run at 

(p)orb ± 2(T. 
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Figure 12. A comparison of the 2002-2003 and 2003-2004 multi- 
site photometric campaigns (MSC) SPB terms (upper panel) with 
the 13 strongest terms from Table 1 having fj < 1.3 d — 1 (lower 
panel). The MSC /j (j = 1, ..,6) terms are labeled in the upper 
panel, and the MOST fj (j = 1, ..,7) terms, in the lower panel; 
the six fainter MOST terms are marked with solid triangles. The 
horizontal line beginning at /j in the upper panel has a length 
equal to 0.083 d — 1 , the frequency resolution of the MOST data. 



traduction), are compared in Fig. [12] with the 13 strongest 
terms from Table 1 having fj < 1.3 d -1 . In the upper panel 
are plotted the i;-filter amplitudes of the /j (j = 1, .., 6) 
terms from table 7 of JHS05, i.e., those obtained from 
the combined 2002-2003 and 2003-2004 MSC data. In the 
lower panel are plotted the MOST amplitudes of the fj 
(j = 1,..,11, 13, and 16) terms from Table 1; the seven 
strongest ones are labeled, the remaining ones are marked 
with solid triangles. The horizontal line beginning at f{ in 
the upper panel has a length equal to the frequency resolu- 
tion of the MOST data. Clearly, the fe, fi, fs, f-2 quadruplet 
would not be resolved in the frequency analysis in Sect. 3.1. 

In the time domain, the situation is the following. The 
/<3 i/ii/s 1/2 quadruplet gives rise to a multiperiodic beat- 
phenomenon. The shortest beat-period is equal to 7.5 d, 
the three longest ones range from 20.9 d to 23.8 d, and the 
dominant beat-period, due to the interference between the 
largest-amplitude terms f{ and f' 2 , is equal to 11.8 d. If 
the quadruplet were present in the light variation of fi Eri 
when the star was observed by MOST, the 12-day light-curve 
seen in Fig. [1] would include one dominant beat-cycle and 
only halves of the longest ones. The question now is whether 
the short time-span of the MOST data or, equivalently, its 
limited frequency resolution, could be the sole reason for the 
appearance of the two frequencies, /1 and fo, in the lower 
panel of Fig. [T^] instead of the quadruplet. 

In order to answer this question we carried out fre- 
quency analysis of two sets of synthetic time-series. Both 
sets consisted of the magnitudes, rrii, computed from the 
equation 

6 

m != ^4 sra(27r/j*< + $;•), (5) 
3=1 

where ti are the epochs of the MOST observations. In set 
1, the amplitudes, Aj, and the frequencies, /j, were taken 



from table 7 of JHS05, i.e., they were equal to those plotted 
in the upper panel of Fig. 1121 The initial phases, , were 
adjusted so that the rms difference between the MOST out- 
of-eclipse detrended light-curve and the synthetic light-curve 
was a minimum. In this way we have shifted the synthetic 
beat-pattern in time to best represent the observed varia- 
tion. Note that because of the lapse of time, the MSC initial 
phases are useless in this context. 

In set 2, only the frequencies, /j, were from table 7 
of JHS05, while the amplitudes, A'j, and the initial phases, 
$<, were adjusted to minimize the rms difference between 
the MOST and synthetic light-curves. In other words, the 
two sets represent error-free 2002-2004 light-curves of \i Eri 
sampled at the epochs of the MOST observations, with ad- 
justed initial phases in set 1, and the amplitudes and phases 
in set 2. 

The results of a frequency analysis of sets 1 and 2, car- 
ried out in the same way as the frequency analysis in Sect. 
3.1, are shown in the upper and lower panel of Fig. 1131 re- 
spectively. In both cases, the seven strongest terms are dis- 
played. As can be seen from the figure, only four terms (solid 
lines) have survived with sizable amplitude. The remain- 
ing three terms (triangles) have amplitudes ranging from 
0.192 mmag to 0.105 mmag in the upper panel, and from 
0.524 mmag to 0.144 mmag in the lower panel. In both pan- 
els, the /g , f[ , f' 5 , f% quadruplet is replaced by two terms 
with the amplitudes and frequencies rather close to those of 
the MOST fi and fi terms. In the upper panel, the first 
term (the highest-amplitude solid line) has an amplitude 
very nearly equal to that of the f\ term and the frequency 
0.013 d _1 greater than f\. The second term (the second 
highest-amplitude solid line) has an amplitude 0.81 mmag 
smaller than that of the fi term and the frequency 0.007 d _1 
greater than f%. In the lower panel, the first term matches 
the /1 term almost perfectly: it has the same frequency and 
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Figure 13. The results of a frequency analysis of synthetic time- 
series having the same frequency resolution as the MOST data 
(solid lines and triangles). The MSC f. (j = 1, .., 6) terms (dashed 
lines) are shown for comparison. Upper panel: the synthetic time- 
series was computed with the MSC amplitudes and frequencies 
but the initial phases adjusted to minimize the rms difference 
between the MOST light-curve and the synthetic one (referred to 
as set 1 in the text). Lower panel: the same for set 2 in which not 
only the initial phases but also the amplitudes were adjusted. In 
the upper panel the dashed lines show the MSC amplitudes (the 
same as the solid lines in Fig. I12[l . while in the lower panel are 
shown the adjusted amplitudes. 



amplitude only 0.26 mmag greater. The second term fares 
only slightly less well: its frequency and amplitude differ 
from those of the fi term by 0.008 d _1 and 0.32 mmag, 
respectively. Thus, the question posed at the beginning of 
this section can be answered in the affirmative: the limited 
frequency resolution of the MOST data is the sole reason for 
the appearance of the two frequencies, /i and fi, instead of 
the fg,fi,fs, fi quadruplet. We conclude that the four fre- 
quencies have not changed appreciably during the six years 
between the epochs of the MSC and MOST observations. 
In addition, the dashed lines in the two panels of Fig. [T3] 
show no amplitude change of f' 5 and f%, but a decrease in 
the amplitude of f{, and an increase in that of /g, indicat- 
ing long-term amplitude variability of these two terms of the 
quadruplet. 

^From Fig. 1131 one can also see that terms f' 3 and f± 
are retrieved almost unchanged. The MOST terms closest in 
frequency to and f'% are ft and /4, respectively. As can be 
seen from Fig. 1121 terms f' A and / 4 are close to each other in 
amplitude and frequency. We conclude that they are almost 
certainly due to the same underlying variation, unchanged 
since 2002-2004. The case of fg and ft is less certain because 
the frequency difference amounts now to 0.022 d" 1 , almost 
twice the difference between / 4 and /J. However, given the 
limited frequency resolution of the MOST data, the possi- 
bility that fs and ft are related cannot be excluded. 

Let us now consider the small-amplitude synthetic 
terms, shown with triangles in Fig. 1131 The highest- 
amplitude term, at 0.525 d _1 in the lower panel of the figure, 
is virtually identical with term /13 of Table 1. Another one, 



at 0.748 d" 1 in the upper panel, is close in frequency to term 
/16 but has an amplitude smaller by a factor of about 0.4. 
The remaining ones bear no relation to the terms listed in 
Table 1. 

9.2 SPB terms added by MOST 

The low- frequency terms fa, fa, fa (solid lines), and fa (the 
triangle at 0.109 d^ 1 ), plotted in the bottom panel of Fig. 1121 
have no counterparts in the MSC frequency spectrum of /1 
Eri. This does not necessarily mean that these terms were 
absent in the star's light-variation in 2002-2004. They all 
have amplitudes smaller than 2.9 mmag, and therefore may 
have been missed by JHS05 because of the 1 // increase of 
periodogram noise in ground-based photometry. 

The three terms seen close to /4 in the lower panel of 
Fig. [E2 viz., /10 at 1.035 d~\ f n at 1.115 d~\ and / 8 
at 1.254 d _1 , have amplitudes below the detection thresh- 
old of MSC. With /4, they form a quadruplet consisting 
of a barely-resolved equidistant triplet fio, fn, fi, and the 
rightmost term fa, separated from fa by about 3/4 of the 
frequency resolution of the data. 

As can be seen from Table 1, all terms with frequen- 
cies higher than 1.3 d _1 have amplitudes smaller than 0.71 
mmag. The two frequencies of largest amplitude are /12 at 
1.825 d" 1 and f lt at 2.423 d" 1 ; /12, /15 at 1.620 d _1 , and 
/19 at 1.722 d" 1 form an equidistant triplet with a separa- 
tion of 0.103 d" 1 , while /14 is preceded by terms /rg and 
/24, and followed by fat and /20. The frequency separation 
between the preceding terms is very nearly equal to that 
between the following terms, and is equal to the separation 
between the members of the /12, /rs, /19 triplet. 

The remaining terms with frequencies lower than 
15 d _1 , viz., / 22 at 3.780 d~\ fat at 4.850 d _ \ and / 3a 
at 5. 499 d~ , have S/N ^ 3.5, and therefore are insignifi- 
cant (|f3reger et al.|[l993T ). Thus, the lowest SPB frequency in 
the light-variation of jj, Eri we have found in the MOST data 
is f 9 = 0.109 d _1 , while the highest one is f 20 = 2.786 d _1 . 
However, keeping in mind the synthetic time-series lesson of 
the preceding section, we warn that neither these frequen- 
cies nor the frequencies in the multiplets around / 4 , /12, and 
/14 should be taken at their face values. 



10 SUMMARY AND SUGGESTIONS FOR 
IMPROVEMENTS 

A frequency analysis of the MOST time-series photometry 
of \i Eri (Sect. 2) resulted in decomposing the star's light- 
variation into a number of sinusoidal terms (Table 1) and 
the eclipse light-curve (Fig. [SJ). New radial- velocity obser- 
vations (Sect. 4.1), some obtained contemporaneously with 
the MOST photometry, were used to compute parameters 
of a spectroscopic orbit, including the semi-amplitude of the 
primary's radial- velocity variation, K\, and the eccentricity, 
e (Table 3). Frequency analysis of the radial- velocity residu- 
als from the spectroscopic orbital solution failed to uncover 
intrinsic variations with amplitudes greater than 2 kms - 
(Sect. 4.3). The eclipse light-curve was solved by means of 
the computer program EBOP (Sect. 5). The solution con- 
firms that the eclipse is a total transit, yielding the inclina- 
tion of the orbit, i, and the relative radius of the primary, 
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r p (Table 4). A Rossiter- McLaughlin anomaly, implying the 
primary's direction of rotation concordant with the direc- 
tion of orbital motion and V TO t sini of 170 ±45 kms -1 , was 
detected in observations covering ingress (Sect. 6). 

^Frorn archival UBV and uvby/3 indices and the re- 
vised Hipparcos parallax, we derived the primary's effective 
temperature, T e g = 15 670 ±470 K, surface gravity, logg = 
3.50 ±0.15, bolometric correction, BC = -1.40 ±0.08 mag, 
and the luminosity, logL/L = 3.280 ±0.040 (Sect. 7.1). 
These values of T e « and logL/LQ place the star in the HR 
diagram f Fig. Ill[> in very near ly the same position as in 
figure 5 of lHandler et all (|2004f ). Thus, these author's con- 
clusion that u Eri is just before or shortly after the end of 
the main-sequence stage of evolution remains unaltered. A 
discrepancy between the effective temperature and bolomet- 
ric co rrection we derived and those obtained from iFlower'sl 
calibration of the B — V index was noted. An analysis 
of a high signal-to-noise ratio spectrogram yielded T e g and 
log<? in good agreement with the photometric values. The 
abundance of He, C, N, O, Ne, Mg, Al, Si, P, S, CI, and 
Fe were obtained from the same spectrum by means of the 
spectrum synthesis method. 

A mean density, (p) ph = 0.0418 ± 0.0057 gem -3 , was 
computed from the star's HR coordinates and the mass 
M p = 6.2 ±0.2 Mq, estimated from an X = 0.7, Z = 0.02 
evolutionary track under assumption that the star is still in 
the main-sequence stage of its evolution. Another value of 
the mean density, (p) or b ~ 0.035 ± 0.003 gem -3 , very nearly 
independent of M p , was computed from Ki, e, i and r p 
(Sect. 8). With the help of evolutionary tracks, this allowed 
placing the star in the HR diagram on a line (p) or b ~ const. 
For evolutionary tracks with X — 0.7, Z = 0.02, this 
line runs close to the star's position mentioned above. For 
M p > 5 Mq , the line crosses the tracks at an advanced stage 
of hydrogen burning in the core, a phase of secondary con- 
traction, or a phase of shell hydro gen-burning (Fi g. mi) . This 
is the same conclusion as that of lHandler et al.l (|2004T ) but 
independent of photometric calibrations and the Hipparcos 
parallax. 

The difference between the two values of mean density 
of the primary component of u Eri, A(p) — (p) p h — (p)orb, 
is equal to 0.0070 ± 0.0065 gem" 3 (Sect 8). The differ- 
ence, although large, is not significant. Thus, the discrep- 
ancy mentioned in the Introduction ha s vanished. It is now 
clear that the value of (p) or b derived bv lJerzvkiewicd {2009), 
viz., 0.0132 ± 0.0050 gem -3 , was too small by a factor of 
2.6. The sole reason for this is that the relative radius of the 
primary component obtained from the ground-based pho- 
tometry, r p = 0.287, is 1.4 times greater than r p from the 
MOST photometry (Table 4). This, in turn, was caused by 
insufficient accuracy of the ground-based photometry. 

The standard deviation of A(p) cannot be easily re- 
duced because it is determined mainly by the standard devi- 
ation of (p) p h which, in turn, is dominated by the standard 
deviation of the luminosity (Sect. 8). Reducing the latter 
would require more accurate bolometric corrections (Sect. 
7.1) or, in other words, better absolute bolo metric fluxes 
of hot stars than those provided years ago by ICode et alj 
1 19761) . As far as we are aware, such data are not forthcom- 
ing. On the other hand, (p) or b can be improved readily by 
obtaining a more accurate value of Ki than that given in Ta- 
ble 3. A radial-velocity project undertaken to achieve this 



purpose should include measurements covering the eclipse 
so that the velocity of rotation of the primary component 
could be better constrained than in Sect. 6. 

The frequency resolution of the MOST observations of 
u Eri is insufficient to resolve the closely-spaced SPB fre- 
quency quadruplet, 0.568 d -1 , 0.616 d -1 , 0.659 d -1 , and 
0.701 d -1 , discovered from the ground by JHS05. Despite 
this, we showed that the four frequencies have not changed 
appreciably during the six years between the epochs of the 
ground-based and MOST observations (Sect. 9.1). The other 
two SPB frequencies seen from the ground, 0.813 d -1 and 
1.206 d - , may also have their MOST counterparts. More- 
over, our analysis of the MOST data added 15 SPB terms 
with frequencies ranging from 0.109 d -1 to 2.786 d ; the 
amplitudes range from 2.843 ±0.013 mmag at 0.387 d -1 to 
0.164 ±0.008 mmag at 2.786 d -1 (Sect. 9.2). However, in 
view of the limited frequency resolution of the MOST data, 
these frequencies should be treated with caution. Finally, 
short-period variations of uncertain origin were found to be 
present in the data (the Appendix). 
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Figure Al. The standard deviations of the residuals from a least- 
squares fit of equation (2) with N = 20 and the parameters from 
Table 1 to the out-of-eclipse data, computed in the data seg- 
ments and the 0.0433-d adjacent intervals of the non-segmented 
data (the left-hand and the right-hand side of the vertical line, 
respectively). The circles denote segments with ~111 data-points 
per hour, and the dots, the segments and the intervals with ~55 
data-points per hour. The pluses indicate the segments and the 
intervals shown in Fig. lA2l 
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APPENDIX A: SHORT-TERM VARIATIONS IN 
THE MOST OBSERVATIONS OF ERI 

As we mentioned in Sect. 2, in many data segments and in 
the non-segmented data there occur short-term variations 
with ranges up to 5 mmag. We shall discuss these variations 
in some detail presently using the residuals, O — C, from 
a least-squares fit of equation (2) with N = 20 and the 
parameters from Table 1 to the out-of-eclipse data. This is 
justified because all terms up to /20 have periods longer than 
0.36 d, about 8 times the length of the longest segment. 
In order to include the non-segmented data, we cut them 
into adjacent intervals of constant length, equal to 0.0433 d, 
the mean length of the segments. Then, we computed the 
standard deviations of the residuals in the segments and the 
intervals. The results are plotted in Fig. I All 

Fig. I A2I shows the segments and the 0.0433-d intervals 



16 



M. Jerzykiewicz et al. 




i 'i 1 1 

141.57 141.58 141.59 141.60 141.61 141.62 148.31 148.32 148.33 148.34 148.35 



HJD-2455000 

Figure A2. The O-C residuals from the N = 2Q fit for the three 
segments (the left-hand panels) and three 0.0433-d intervals (the 
right-hand panels) indicated with pluses in Fig. I All Plotted are 
the data (dots) and the N = 31 synthetic light-curve (lines). 

indicated with pluses in Fig. I All We selected the segments 
and the intervals having the largest, median, and smallest 
standard deviations (from top to bottom, respectively). 

In Sect. 3.1, we mentioned that formally significant 
terms (i.e., those having S/N > 4.0) with fj > 15 d -1 
have frequencies which are either close to whole multiples 
of the orbital period of the satellite or differ from the 
whole multiples by ~1 d _1 . Indeed, denoting the orbital 
period of the satellite by f we have /21 « 5/„+l d _1 , 
/as « 2/ , fa « f +l d-\ fa « 4/ +l d" 1 , / 28 w 3/ D , 
/29 ~ 4/ -l d _1 , and / 30 « 6/ + l d" 1 . Although it is 
not clear why a frequency of ~1 d" 1 should appear in the 
satellite data, all these frequencies are probably instrumen- 
tal. The remaining j > 20 terms have S/N < 4.0, but were 
included in order to remove as much out-of-the-eclipse vari- 
ation as possible from the light-curve. From Fig. IA2I it is 
clear, however, that the N — 31 solution leaves a part of 
the short-period variations unaccounted for. Adding further 
periodic terms did not help. Thus, the unaccounted-for part 
must consist of incoherent short-period variations. Whether 
they are intrinsic or instrumental is impossible to say. 

This paper has been typeset from a TgX/ I^TgX file prepared 
by the author. 



